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ABSTRACT 

Fermi acceleration can take place at ultra-relativistic shock waves if the upstream or down- 
stream magnetic field has been remodeled so that most of the magnetic power lies on short 
spatial scales. The relevant conditions under which Fermi acceleration become efficient in the 
presence of both a coherent and a short scale turbulent magnetic field are addressed. Within 
the MHD approximation, this paper then studies the amplification of a pre-existing magnetic 
field through the streaming of cosmic rays upstream of a relativistic shock wave. The mag- 
netic field is assumed to be perpendicular in the shock front frame, as generally expected in 
the limit of large shock Lorentz factor In the MHD regime, compressive instabilities seeded 
by the net cosmic-ray charge in the shock precursor (as seen in the shock front frame) develop 
on the shortest spatial scales but saturate at a moderate level dBjB ~ 1, which is not sufficient 
for Fermi acceleration. As we argue, it is possible that other instabilities outside the MHD 
range provide enough amplification to allow successful Fermi acceleration. 
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The physics of collisionless shock waves is of paramount impor- 
tance in modem astrophysics, as it governs a variety of phenomena, 
most notably the emission of high energy radiation. In particular, 
the afterglow radiation of gamma-ray bursts is generally attributed 
to the synchrotron emission of electrons that have been accelerated 
around the relativistic external shock wave, with shock Lorentz fac- 
tor r^h ~ 100 (Paczynski & Rhoads 1993, Katz 1994, Meszaros & 
Rees 1997, Sari & Piran 1997, Vietri 1997, Waxman 1997; see Pi- 
ran 2005 for a review). In this framework, it can be shown that 
the magnetic field, downstream of the shock wave, must have been 
amplified by orders of magnitude beyond the shock compression 
of a typical interstellar magnetic field (see for instance Waxman 
1997; see also Piran 2005 and references therein). Recent work by 
Li & Waxman (2006) further shows that even the upstream mag- 
netic field must have been amplified by at least one or two orders 
of magnitude in order to explain observed X-ray afterglows. 

These results should be put in perspective with recent studies 
on the process of Fermi acceleration in the test particle limit around 
relativistic shock waves. In particular, Niemiec & Ostrowski (2006) 
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have shown that Fermi acceleration becomes inoperative in a turbu- 
lent magnetic field with power spectra of the Kolmogorov or scale 
invariant type. Through analytical calculations, Lemoine, Pelletier 
& Revenu (2006) have shown that relativistic Fermi acceleration 
is inefficient if the magnetic field power is distributed on scales 
larger than the typical Larmor radius of the accelerated population, 
in agreement with the above numerical result. Therefore the inter- 
pretation of the gamma-ray burst afterglow emission from the rel- 
ativistic external shock wave requires the magnetic field to have 
been amplified on short spatial scales. 

It has been suggested that the relativistic two stream Weibel 
instability could amplify the downstream magnetic field to the level 
required by afterglow modeling of gamma-ray bursts (Gruzinov 
& Waxman 1999; Medvedev & Loeb 1999). Subsequent studies 
have however argued that this instability should saturate early on 
(Wiersma & Achterberg 2004; Lyubarsky & Eichler 2006). The 
question of the long term evolution of the downstream magnetic 
field (on timescales » w^' , with cOp the plasma frequency) also re- 
mains open. Ongoing particle-in-cell simulations should eventually 
shed light on this issue (Silva et al. 2003; Frederiksen et al. 2004; 
Medvedev et al. 2005; Kato 2007; Chang, Spitkovsky & Arons 
2008; Keshet et al. 2008; Spitkovsky 2008). 

Regarding the amplification of the upstream magnetic field, 
the present situation is reminiscent of results obtained for non- 
relativistic supernovae remnant shock waves, where the interstellar 
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magnetic field has apparently been amplified by one or two orders 
of magnitude (see Volk, Berezhko & Ksenofontov 2005, Parizot 
et al. 2006). In this case, the leading candidate for the instability 
is the streaming instability, seeded by the cosmic ray precursor in 
the upstream plasma (see Bell 2004, 2005; Pelletier, Lemoine & 
Marcowith 2006, Marcowith, Lemoine & Pelletier 2006; Reville et 
al. 2007; Niemiec et al. 2008; Amato & Blasi 2008; Zirakashvili, 
Ptuskin & Volk 2008; Reville et al. 2008). The generalization of 
this instability to the relativistic regime has been studied on phe- 
nomenological grounds by Milosavljevic & Nakar (2006) who have 
concluded that it should be able to account for the degree of ampli- 
fication inferred in gamma-ray bursts. More recently, Reville, Kirk 
& Duffy (2006) have derived in detail the dispersion relation for a 
parallel shock wave and the saturation due to thermal effects. 

In the present work, we propose to explore in detail the gen- 
eralization of this type of instability to the ultra-relativistic regime 
Fsh » 1. One crucial difference with the previous works on the 
streaming instability is that we consider the most natural case of su- 
perluminal shock waves (with a magnetic field perpendicular to the 
shock normal in the shock front frame). This case is more generic 
than the parallel configuration studied previously because the trans- 
verse component of the magnetic field is boosted by the shock 
Lorentz factor when going to the shock frame. Another important 
difference is that we bring to light a new type of instability, of a 
compressive nature. Finally, in contrast with most particle-in-cell 
simulations performed to date, our study focusses on magnetized 
shock waves, for which there exists a coherent upstream magnetic 
field (whose dynamical influence on the shock jump conditions 
can be neglected however). Nevertheless there exist pioneering PIC 
simulations in the moderately relativistic regime, which include 
both a mean field and a significant mass ratio between electrons 
and ions, see Hededal & Nishikawa (2005), Dieckmann, Shukla, & 
Drury (2008). 

We adopt a simplified description in which the cosmic-ray dis- 
tribution is modeled as a step function out to some distance („ and 
we neglect the cosmic-ray response to the disturbance. This latter 
assumption is justified by the fact that the instability is maximal on 
the shortest spatial scales, orders of magnitude below the typical 
Larmor radius of accelerated particles. 

The paper is organized as follows. In Section[2] we introduce 
the main scales of the problem, most notably the diffusion scale of 
the cosmic rays; we then calculate the level of amplification that 
is necessary to make Fermi acceleration operative. Section|3]is de- 
voted to the investigation of the instabilities under the condition 
that some cosmic-rays have undergone a first Fermi cycle. Section|4] 
summarizes out results and provides some outlook. Details of the 
calculations are provided in Appendix IaI 



2 GENERAL CONSIDERATIONS 

We carry out most of the discussion in the shock front rest frame, 
hence unless otherwise noted, all quantities are evaluated in this 
frame. We use the subscripts in or |d to tag quantities measured in 
the upstream or in the downstream rest frame respectively. 



2.1 Upstream diffusion lengtfi 

In the upstream rest frame, cosmic rays can never stream too far 
ahead of a relativistic shock wave since this latter propagates to- 
wards upstream with velocity Vji, = p^^^c « c [the shock Lorentz 
factor Fsh = (1 "Z^sii) ''^ ^ ^l- Cosmic rays scatter on magnetic 



turbulence upstream before they are caught back by the shock wave 
when their pitch angle Sju ~ 1/Fsh (Gallant & Achterberg 1999, 
Achterberg et al. 2001). Consequently, they can travel a distance 
4t|u, which may take the following values depending on the ratio 
of the Larmor radius tliu to the coherence length of the upstream 
magnetic field A^^^ (Milosavljevic & Nakar 2006): 



for small scale turbulence 



''Criu 






for large scale turbulence 



(?'L|u S> Fsh/ic|u) 



(''Liu ^ Fsh/ic|u) 



(1) 



(2) 



Both regimes, short or large scale turbulence can be expected 
at some point, insofar as the excitation of the upstream magnetic 
field on short spatial scales is due to the streaming of the non- 
thermal particle population in the shock precursor. Indeed, cosmic 
rays of the first generation are to interact with a turbulent magnetic 
field ordered on large scales. However, provided the instability that 
they trigger grows fast enough, cosmic rays of the next generation 
will propagate in short scale turbulence. In reality, the situation is 
likely to be more complex as the process of particle propagation up- 
stream and magnetic field generation are closely intertwined. The 
fact that the non-thermal population contains particles of different 
energies, which can stream at different distances from the shock 
front, should also play a significant role. In this respect, Keshet et 
al. (2008) have observed that the upstream magnetic field is af- 
fected to greater distances as time goes on. This strongly suggests 
that higher energy cosmic rays are produced as time goes on, and 
that, by travelling farther in the upstream medium, they excite the 
turbulence at larger distances from the shock front. 

In the discussion that follows, we estimate the growth of un- 
stable modes in both the limits of small or large scale turbulence in 
order to remain as general as possible. Certainly the limit of large 
scale turbulence is more restrictive with respect to the growth of 
the instability, since the distance travelled upstream is significantly 
reduced with respect to that in small scale turbulence. Whichever 
limit prevails depends on the ratio of the short scale turbulent to the 
large scale (coherent) magnetic field strength, as discussed in Sec- 
tion |2.2| We also discuss the effect of higher energy cosmic rays on 
the growth rate. 

It is important to emphasize that the distance that controls 
the growth of the instability is that between the shock front and 
the position of the particle, which is smaller than 4i|u by a factor 

(1 -fish) ~ (2Fs|,) • In the following, we will need the expres- 
sion for 4r|sh (also noted 4r)i i-e. the length scale of the cosmic-ray 
distribution as measured in the shock front rest frame. It can be cal- 
culated by transforming the upstream residence time /.lu ^ 4i|u/f 
in the shock front frame /,.|sh = ?,|u/Fsi,, and then by rewriting in 
the expression obtained the upstream Larmor radius and coherence 
length in terms of their shock frame equivalent. In the perpendicu- 
lar (or superluminal and ultra-relativistic) configuration of interest, 
''Lish - F^ii^rLiu, /lc|sh - F^|j'/ic|u. This boost of the coherence length 
is valid for wavenumber modes that are parallel to the shock nor- 
mal; perpendicular wavenumber modes remain unchanged. In the 
following, we thus use the following expressions for I^y'. 



for small scale turbulence 



f = f 



cr|sh 



'L|sh 



(3) 
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for large scale turbulence 



4r 



aish 



(4) 



Below, we find a process of generation of intense magnetic 
disturbances at short spatial scales. The expected instabilities in all 
relevant cases generate disturbances with a coherence length Ac\a 
between the minimal length /mhd = c/aipj required for the validity 
of MHD description and the diffusion length 4r|ii of the cosmic 



rays, but with a preference for short scale, i.e. /mhd ^ 'Iciu 
The minimal MHD length is: 



< e, 



'mhd = p. 



A'Olii ■ 



cr|u* 



(5) 



where /?a = So|u/ ^/4jrpu\uC- denotes throughout this paper the 
Alfven velocity (measured upstream) in units of c, and ro|u = 
mpC^/eBoiu is the Larmor radius of thermal protons in the upstream 
comoving frame. The above MHD scale is measured in the up- 
stream plasma rest frame. 



2.2 Requirements for successful Fermi acceleration 

As discussed in Lemoine, Pelletier & Revenu (2006), a particle can 
execute a large number of Fermi cycles through the shock only if 
the turbulence coherence scale is much smaller than the Larmor ra- 
dius (assuming no coherent magnetic field), in detail r^^ s> /iciursh 
upstream, or r^^j » /l^d as measured downstream. 

The Fermi process becomes inoperative in the limit of large 
coherence length (when compared to the typical Larmor radius) 
since the cycle time becomes smaller than a coherence time hence 
the field is effectively regular and superluminal on a typical Fermi 
cycle. Conversely, the Fermi process will be operative provided ei- 
ther of the above inequalities is fulfilled, since the memory of the 
initial conditions at shock crossing is then erased by pitch angle 
diffusion in the short scale turbulence. 

In the generic case of a ultra-relativistic superluminal shock, 
the above two conditions are actually related by the shock jump 
conditions, since the Larmor radius increases by ~ F^^j when go- 
ing from downstream to upstream, while the coherence length in- 
creases by ~ Fiei. The quantity Fiei corresponds to the Lorentz fac- 
tor of the downstream fluid as measured upstream, and in the case 
of a ultra-relativistic strong shock, Fjei - Fsh/ yl (Blandford & 
McKee 1976). 

In the following, we describe the upstream and downstream 
fluids as comprising a large scale component Bo and a short scale 
turbulent generated by some unspecified instability. In this case, 
Fermi acceleration will be efficient provided the short scale field 
has a sufficient amplitude as compared to Bq. We discuss this con- 
dition in the following. 



2.2.1 Upstream motion 

In this particular section, all quantities are evaluated in the upstream 
rest frame and we avoid using the notation [„ out of clarity. We de- 
rive here the requirements on the amplitude of the short scale com- 
ponent of the magnetic field that would allow successful Fermi ac- 
celeration. Without loss of generality, we assume that the coherent 
component Bo lies in the plane x-y, the x direction corresponding 
to the shock normal (oriented toward upstream infinity). We de- 
note by b the irregular component (whose amplitude is expressed 
in units of 5o)- The layout is pictured in Fig.[T] 

The total field is thus written as B = Bo(e + h), and the particle 



downstream 




^ y^ upstream 



Figure 1. Schematic plot of the geometrical configuration in the upstream 
plasma rest frame. The coherent component Bo lies in the x - y plane, with 
,v pointing towards the shock normal, and (9, 0) the angles of the particle 
velocity in this frame. 



trajectory is governed by the following equation: 

dJS 

— = WL.o;8x(e-i-b), 
at 



(6) 



with wl.o denoting the Larmor frequency expressed relatively to 
6o- We neglect the influence of any short scale time varying electric 
field tfE in the above equation of motion since (5ii/(So^) ~ w/fc <c 
1. 

In order to quantify the return timescale and the condition for 
successful Fermi acceleration as a function of the amplitude of the 
random component of the magnetic field, it is useful to express 
this equation in terms of the angles Q and (f), which are defined as: 
y8_v = cosd,/3y = sin 6 cos and/?. = sin 6 sin <^. These equations of 
motion read: 



9 = aiL.o \(ey + by] sin cp- b- cos (f>\ , 

/ x cos d> 
- fc, + e, + by + b. 



wl.0 



cos<p . smi^ 
tan0 "~tane 



(7) 
(8) 



Rescaling the Larmor frequency with respect to eyBo, i.e. defining 
OIL ,, = eyCiJi^Q, and noting that <c 1 when the particle propagates 
upstream, the equations of motion for 9 and can be rewritten as 
follows: 



coi^y sin(0) -I- (i>L,yA sin(0 - o") , 
cos cp cos(9 - o") 



9 



9 



(9) 



The factors A and cr describe respectively the amplitude and the 
phase of the small scale magnetic field in the shock front plane, i.e. 
by = A cos cr and b-- A sin a (assuming isotropic turbulence). The 
amplitude A = SB/ By is measured with respect to By. 

One can extract from the above system the unperturbed trajec- 
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tory, i.e. assuming A = 0: 



COS0 = COS(l>i 



sin 6*1 



(10) 



sin 6 

where values indexed with -, are calculated at some initial time. 
Since entry into upstream corresponds to cosSj > ySsh hence 
Oi ^ l/Fsi,, and exit from upstream corresponds to 9 > l/Fsi,, the 
above equation suggests that <p is driven to an angle close to ±n/2. 
The sign is given by the sign of e,,: for e,, > 0, it is +n/2. This 
region lies opposite to that which allows return to the shock when 
the particle travels downstream (see Lemoine, Pelletier & Revenu 
2006), hence Fermi cycles cannot be completed. It is important 
to note that this unperturbed trajectory is executed on a timescale 
funpert - '■L.,./(FshC) (where rL,v = c/aiL,,). 

The noise term comes from the random phase cr (assuming 
isotropic turbulence in the shock front plane, an assumption that 
can be relaxed). Therefore the trajectory deviates from the unper- 
turbed trajectory to linear order by (Aip, AG), whose variance in- 
creases linearly with time: 

TcAr , 



<A<^2) 



3^ 



02 



{AO^} 



■ A^WL.vT^Af . 



(11) 



Note that the above equation does not contain all the terms 
driving the variance of (p and 9: it neglects term of the form 
f dtidt2{A(p(ti)A(p(t2)} and similarly in A9. However, these terms 
are smaller by 0(A^) as compared to those above, hence they can 
be neglected in a first approximation (we will find A » 1 further 
below). 

In the absence of the regular driving term in the equations 
of motion (i.e. when By = 0), return occurs over a timescale 

fdiff - (rsh'^^'^L •^c) • Noise dominates over the unperturbed tra- 
jectory if (A^^) > 9^, where 9 indicates the unperturbed trajec- 
tory. This inequality must be satisfied within a (return) timescale 
funpert Otherwise the particle will have exited from upstream (along 
the weakly perturbed unperturbed trajectory) before the noise could 
overcome the unperturbed driving force. This implies: 

A^ > . (12) 

FshWLjTc 

If one defines the Larmor radius cl measured with respect to 

the total magnetic field, tl = c\{A^ + l)''2aiL,o , hence tl - 
c/(A(i;L,y) when A » 1, the former inequality amounts to: 

rL 



A > 



Fsh'lc 



(13) 



If this inequality is satisfied, then one can check that (A^^)''^ ~ 
(9(1) on a return timescale fdiff, which implies that the return di- 
rections are isotropized in the shock front plane. Fermi accelera- 
tion should then be efficient (provided the return probability to the 
shock as calculated downstream is also isotropic in (p, see below). 

An interesting implication of the above is to restrict Fermi ac- 
celeration to a rather limited range of Larmor radii: 



Tsh'tc < f^L < AF^h-lc 



(14) 



Very high values of the amplification factor are thus required to 
produce powerlaw spectra over a large dynamic range. 

2.2.2 Downstream motion 

All quantities in this subsection are evaluated in the downstream 
rest frame. In this frame, the coherent magnetic field is nearly fully 



aligned with the y-axis, and we assume indeed that Bo = Soy. 
It is useful to study the evolution of the velocity components per- 
pendicular (^x) and parallel (jSy) to the direction of Bq: P\\ = Py, 
fi\ = p\+ 01. The equations of motion for the velocity p can then 
be written: 

/3x = <^L.o^ (P.b,-Pzb.) 
P± 

/3|l = couo(Pzb,-p,h,) . 



(15) 



One can check that/3n8||+/3xj0± = as it should since /J is conserved. 
The unperturbed trajectory for the couple of variables (ySy, p^} is 
trivial, i.e. /3x,o = const. = p±j and Pio = const. = /Jn;. Return 
may occur or not, depending on these initial velocity components; 
if it occurs, it does so on a timescale funpen - rj^.o/c (Lemoine, 
Pelletier & Revenu 2006). 

Since b^ is not compressed through shock crossing, b^/b, ~ 
l/Fsh for isotropic upstream turbulence, hence this component 
can be neglected. Integrating the above equations to calculate the 
squared perpendicular displacement, one finds: 



(Axi) 



fdf, { dt2{PAh)PAt2)) 

Jo Jo 

w^oji I dfi I dt2 I df3 I d«4 

Px,i Jo Jo Jo Jo 

P..oi.h)P..o{U){b-Xti)b,it^)) 

C C C A^ 

dti dt2 df3— Te 

Jo Jt, Jo 3 



Pi 



Px., Jo 



:PUh) 



(16) 



In the above equation, Px,o(t) represents the unperturbed trajectory: 
P-1.0 = Pxi cos(ajL,oO - Pzi sin(wL,oO- The function P{t) contains 
powers of t up to t^ , as well as sine and cosine functions of 2a»L,of; 
its dimension is that of c^co^- Noise will then dominate over the 
unperturbed trajectory x^fy = Pi_jCt over a unperturbed return 
timescale if: 



A^ > 



I 



WL,0Tc 



or, equivalently: 



A > 






(17) 



(18) 



As before, tl denotes the Larmor radius measured with respect to 
the total magnetic field. This condition on A is very similar to that 
obtained upstream, since tliu ~ rsh''L|<i and /i^u ~ Fsi,/i;.|(], while 
All, ^ A|d. 

According to the above analysis, and following Lemoine, Pel- 
letier & Revenu (2006), phase mixing should be sufficiently large 
to erase all dependence on the phase of the velocity vectors, hence 
Fermi acceleration should be fully operative, if either one of the 
above conditions on A is satisfied. 

It is interesting to discuss the above results in light of recent 
Monte Carlo numerical simulations of relativistic Fermi accelera- 
tion in the presence of short scale turbulence (Niemiec, Ostrowski 
& Pohl 2006). These authors have investigated the efficiency of 
the Fermi process for various turbulence configurations, including 
a coherent magnetic field, a long wave turbulence and a short scale 
component. They confirm that Fermi acceleration is not operative 
if SB IB = 0, but find that a spectrum can develop over about two 
orders of magnitude when SBjB = 80 (see their Fig. 1). The high 
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energy break can be directly interpreted as the energy at which the 
above inequality Eq. jlSt is no longer fulfilled. Beyond that point, 
the spectrum steepends and the Fermi process becomes inopera- 
tive. Interestingly, these authors also show that if the pitch angle 
scattering amplitude in the short scale turbulence is independent 
of energy, the break disappears and the spectrum continues without 
bound. This is also expected, insofar as the break that is determined 
by Eq. l llSt stems from the fact that the scattering time in the short 
scale turbulence scales as the square of the Larmor time whereas 
the unperturbed trajectory in the background field depends linearly 
on the Larmor time. Assuming a constant pitch angle scattering 
amplitude in the notations of Niemiec, Pohl & Ostrowski (2006) 
implies that the two scattering timescales evolve in the same way, 
hence if short scale turbulence suffices to isotropize directions at a 
given energy, it will do so at all energies. 

In order to produce powerlaw spectra over a large dynamic 
range, it is thus necessary to reach an amplification factor 6BIB as 
large as possible, since this ratio determines the dynamic range of 
particle energies. Interestingly, the afterglow modeling of gamma- 
ray bursts suggest that indeed, the usptream and downstream mag- 
netic fields have been amplified by a large factor. In the range of 
energies in which Fermi acceleration is operative, it is likely that 
the spectral index would be equal to the so-called canonical value 
s = 2.3, although the generality of this prediction for different 
shapes of the 3d turbulence spectrum remains to be studied. The 
results of Niemiec, Ostrowski & Pohl (2006) cannot be used to in- 
fer this spectral index, since they have restricted their analysis to 
values of SBjB < 100, and the spectral indices they report have 
been measured beyond the break energy. Techniques developed in 
Lemoine & Pelletier (2003), Lemoine & Revenu (2006) are partic- 
ularly suited to study this problem and calculations are underway. 



3 INSTABILITIES AT PERPENDICULAR SHOCK 
WAVES 

The magnetic field in the shock front B^^y, is related to the magnetic 
field in the upstream frame 6|„ and the associated electric field Ei^ 
through the Lorentz transform: 



6|l|sh 

B±|sh 



6|l|u 

Tsh (Bx|u -Ah X Ex|u) 



(19) 



To zeroth order, the electric field in the upstream plasma frame van- 
ishes, hence the magnetic field in the shock front frame is mostly 
perpendicular unless B[ii is aligned along the shock normal to within 
an angle ~ 0(1/Tsh) (subluminal shock). It then suffices to consider 
the fully perpendicular situation with B^^^sh = 0. 

The cosmic rays stream ahead of the shock wave, carrying 
a net charge density p^r which will induce a counteracting charge 
density Ppi in the background plasma. Note that the cosmic rays do 
not induce an electrical current at zeroth order in the shock front 
frame, only a net charge density. Since we consider the generation 
of short scale turbulence, we neglect the cosmic ray response to 
this short scale turbulence other than its eifect on the cosmic ray 
distribution scale 4r- For simplicity, we approximate the cosmic- 
ray charge profile with a top-hat distribution. 

As usual, we search for a stationary regime which serves as a 
basis for perturbating the equations in the time-dependent regime. 
The details of the calculations are provided in Appendix lAlfor both 
stationary and time dependent quantities. In particular, the station- 
ary regime exhibits the following set-up: 

B = Byty, u = u^ex + u-Ci . (20) 



Recalling that the shock normal is directed along x, Uj ^ -FshySsh 
characterizes the velocity of the upstream inflowing through the 
front. The velocity component u, along the front is small compared 
to Mt, but its shear d^u^ cannot be neglected (see Appendix lAt. 



3.1 Linear analysis of the reduced system 

The complete system that governs the linear evolution is of seventh 
order (see Appendix [Ajl and thus rather involved. Nevertheless, it 
remains possible to derive the main results since the Lorentz trans- 
form from the upstream comoving frame to the shock front frame 
dominates the MHD propagation effects when the wave scale is 
"not too small", as will be made more precise in a next subsection. 
In other words, the instability growth rates can be derived to lowest 
order by assuming vanishing Alfven and sound velocities. We have 
verified that the analysis of the system including /Ja and jS, (but 
neglecting terms in ii-) does not modify the growth rates obtained 
further below. 

One should point out that electromagnetic waves with disper- 
sion relation a;|„(k|u) are Lorentz transformed into waves with dis- 
persion relation: 



I 1 sh 



(21) 



In particular, for Alfven waves, the next term on the r.h.s. is 
PAkyc/Tsh- Note that the Alfven velocity is always expressed in 
the upstream plasma rest frame, i.e. /Ja = ^viu/ ^^Tcp^c^. There- 
fore it suffices in what follows to consider the limit k- — > 0. The 
limit ky — > can also be considered but it appears more restrictive, 
because it limits the analysis to the evolution of upstream magne- 
tosonic modes. We will thus consider both cases ky = and ky 
finite. Of course, one can use Fourier analysis without mode cou- 
pling as long as the spatial dependence of the system coefficients 
can be neglected, in agreement with our previous approximations, 
and in particular, that the cosmic ray charge density is modeled as a 
step function. In the shock front rest frame, the physical picture of 
the instability is then as follows. Impingent electromagnetic waves 
propagate in vacuum beyond the length scale 4r which character- 
izes the spatial distribution of cosmic rays upstream of the shock 
front. For < x < 4r, the presence of the cosmic-ray charge den- 
sity induces a short scale instability. Hence we consider electro- 
magnetic waves and seek a solution in Fourier space with uj set by 
its vacuum dispersion relation, but with a complex k^ characterizing 
the amplification in the charge layer < jc < fcr- 

In this perpendicular configuration, the return current (in the 
upstream frame) is not responsible for a supplementary tension ef- 
fect, but for a supplementary compression effect. A by perturbation 
leads to a supplementary vertical (i.e. along z) compression that can 
push like the kinetic compression. If we consider a spatial modula- 
tion along the mean field {ky # 0), a vertical perturbation b, gener- 
ates a supplementary compression in the direction of the mean field 
that can push in phase with the kinetic compression as well. 

For large Fsi, even with k^ i^ 0, the system reduces to fourth 
order and can be expressed as the coupling between the propa- 
gation of the vertical perturbed motion 6u^ and the propagation 
of kinetic compression Swiw, where w = p^c^ is the relativistic 
proper enthalpy density of the cold upstream plasma. The system 
of the two coupled equations can be expressed using the system of 
Eqs. l |AT8t ,( |AT9] l,l [AlB in the hmitySA -^ 0,/3, ^ 0. In particular 
in this limit, one notes that fc, a 0. 

Then one obtains a single equation for dii- (see Ap- 
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pendix lA2.2t where derivatives d, cancel out and where the deriva- 
tive of u- is inserted, its second derivative being neglected (we also 
assume ySsh = 1 in the coefficients): 

E^Su, - {ic + u-Kd,) D^Siu + T,yd,b6u, + Tl^^K^dlSu, = . (22) 

The differential operator D is defined by D = Fsh (c"'(9i -ySshSv). 
For a detailed understanding of the instabilities one has to notice 
that the charge is not completely screened by the plasma as a sup- 
plementary electric field is generated along the flow (i.e. oriented 
towards -x): £,- = u,By/Tsh, along with a sheared vertical motion 
M, such that d^u- = k (see Appendix|All. The instabilities stem from 
this sheared motion. The parameter k, which carries the dimension 
of a wavenumber, is defined in Eq. JA23t . To our present order of 
approximation, it can be expressed as: 



TshW 



(23) 



Note that both p^i and By are here evaluated in the shock front 
frame. In the front frame, the parameter k appears divided by Fsh, 
hence we introduce k, = /f/Fsh. 

The reduced form of the equation (see Appendix|A2J2} clearly 
shows the ordering (by setting D = k^flshTshD, Su- = /?shFshMc): 



D U. - -r + 



k:d, - kl . 
, iD'ii. H — -^Du. + —rdyii. 



: (24) 



The inhomogeneity is contained in k^ that vanishes at large distance 
from the shock. At infinity DSu- = 0, which implies that for a mode 
having a specified k-,, the pulsation a) = p^hk^c, as expected. In the 
precursor where k^ t 0, we follow this mode characterized by its 
frequency and look at the modification of its spatial behavior. As 
explained above, we solve this equation by setting c^'dt i--> i[}^\Jcx, 
dx I-* ikxil - s), so that D i-» ikxPstX^hE (D i-» is), where s is 
a complex number, the imaginary part of which characterizes the 
growth rate. It is obtained by solving the algebraic equation: 



4 kl 2 kl klkj. .u, k.^ . 

E +2— e^ -E -J- +; (1 - e)£- 

/C" K ICl 1 sh '^x 







(25) 



Actually, because e is small, the equation can be simplified to: 

. kl. m 

s^ — '- 

kl 



--r- 



0. 



(26) 



The contribution in u, introduces a small correction to the real part 
of the roots (most u-^ contributions can be canceled out by a frame 
transformation, but that correction cannot be canceled out because 
it corresponds to the electric field component £,). 



3.1.1 Analysis of the case ky = 

Let us first analyze the case ky = 0, or less restrictively for kl <c 
Ifi^l kl. This is a pure magneto-sonic compression with b along the 
mean field. One gets the following three roots: 



-(r-. 



1±/V3 



(27) 



One obtains an instability spatial growth rate y,. = Im(/:,.e) that 
increases with large kx. 



7-v 



V3 



k^J^kl'^ {kx » k.). 



(28) 



Thus the characteristic wave number associated with the charge in 
this problem of magneto-sonic wave propagation is k, = /c/Fsi,. 



3.1.2 Analysis of the case ky ^ 

; recove 
while for ky » k, kj^ one derives : 



With ky i^ 0, one recovers the previous growth rate if ky <k k, kj , 



1/2 



^1 (h-hi,-i).. 



or, equivalently: 



vI/2 



JX - (^,1.)'" {ky » kl'^k^) 



(29) 



(30) 



The spatial growth rates are thus y,. ~ (klkxY^^ and (fc.^,,)''- 
respectively. Remarkably they are independent of k that can be 
chosen arbitrarily, provided it remains small compared to kx for the 
consistency of the derivation. 

3.1.3 Comparison to MHD and cosmic-ray length scales 

One should compare the scale defined by k, to the smallest scale 
'mhd for our MHD description. Note that this length scale is defined 
in the comoving upstream frame. Since kx\u - A;t/Fsi, while ky^^ ^ 
ky, it suffices to require A;./mhd < 1 to ensure that there exist modes 
in the MHD range with k > k,. One obtains: 



k,L 



rlw 

sh 



'3Aro\u 



(31) 



which can be rewritten in terms of the fraction ^^i of (downstream) 
shock internal energy converted into (downstream) cosmic ray en- 
ergy: 



^cr = 



4riw 



(32) 



To this effect, we assume that the accelerated population can be 
described as a power-law of index -s and minimum momentum 
Pmin, so that: 

11 - ii e„ 

"cr - ^ \^^. (33) 

|2 - S\ PminC 

This equation assumes .s > 2; if .s = 2, then |2- .s| should be replaced 
by log(pmax/Pmin). Since e„|d ~ ecrish (as Fshjd = VWS for a strong 
ultra-relativistic shock), one obtains: 

Jl 



k,lf. 



1^ ''I 'min|u 



(34) 



Here r. 



minlu 



denotes the minimum Larmor radius of the accelerated 



population, as measured upstream. The typical energy of the first 
generation of cosmic rays is F^j^mpC', so that r,nmiu ~ r\>'o|u. There- 
fore 



kJu 



1^^ 
\2-s\ 



^crySA « 1 • 



(35) 



Thus k, defines an MHD scale. 

In order for the instability to be efficient, one must also re- 
quire y,4r ^ 1, which is a non-trivial requirement in view of the 
restricted value of t„. One finds: 



7-v4r 



4i;^2^,„,, 



(36) 



|2 - s\ k, 

with g = 1 for large scale turbulence, and g = rg^^/A^^^ for short 
scale turbulence. In the latter case, one may assume A^i„ ^ Imhd, 
as suggested by the fact that the growth rate increases with k, so 
that g ~ l/ySA- Then, in both cases, the total growth is presumably 
much larger than unity if fcr is not too small and kx, ky sufficiently 
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large as compared to k,. Hence even the first generation of cos- 
mic rays is able to destabilize the upstream magnetic field on the 
shortest scales. The effect of the high energy part of the accelerated 
population will be addressed shortly. 

As mentioned, the study of the more extended system includ- 
ing the contribution of the terms involving the Alfvenic and sonic 
contributions (but neglecting terms of order ii^ in the system) con- 
firms the results above for the growth rates, provided /3a ^ 1 and 
A « 1- 



can only be present if acceleration to smaller energies has been 
completed. However it is easy to verify that, if highest energy cos- 
mic rays are present in the shock precursor, they are bound to play 
a dominant role in the amplification of the magnetic field, most 
notably because the length scale of the cosmic ray distribution in- 
creases as r^. For instance, at a distance x from the shock front 
where only cosmic rays of momentum larger than p, can be found, 
the critical wavenumber for the instability reads for s > 2 (s spec- 
tral index of the accelerated population): 



3.1.4 Magneto-sonic saturation 

One can estimate the amplitude of the various components in terms 
of by. It has already been mentioned that 1^,1 ^ \by\. 



k. = k. 



and 



by ^ -^ou^ 

1 sh 



— dySu^ 



Thus we find that b- may achieve a large amplitude, because 



kPa 



dyby . 



(37) 



(38) 



(39) 



In the frame of linear theory, the saturation level of these magneto- 
sonic instabilities is simply estimated by the fact that the com- 
pression has a limited amplitude |(5w/Mi < 1. Defining the power 
spectrum of by and b, per log interval of wavenumber Vt^. = 
(fe^/2;r^)|fcyp in terms of the Fourier component by and similarly 
for Pi,^ in terms of b-, one derives from Eq. l lA33t : 



by 



and 



■P 



1/2 



\E^kl + 2kl\ 






(40) 



(41) 



Note that the power spectrum is normalized to the ratio of ampli- 
tude of the turbulent (small scale) component to the coherent Bo- 
Consider the first branch of the instability, given by Eq. J28I ) which 
applies in the limit k^ » k,, ky <c kj kj . There one can show 
that f ^'- saturates at a value of order {kJk.Y^'^ <& I, but pI''^ 

saturates at a value of order unity when ky ~ kl k'J . For the 
second branch, given in Eq. BOb and which applies in the limit 
ky » k, kj^ , one finds that f^y saturates at a value of order 

k,lky <c 1, but again Pj' saturates at a value unity. This moderate 
saturation level can be directly traced back to the saturation of these 
compressive modes at the linear level |(5w/w| ~ 1. As discussed in 
Sec. 12.21 such a level dB/Bo ~ I should not allow a powerlaw spec- 
trum to develop through Fermi acceleration. 

Nevertheless, our present analytical description is by defini- 
tion limited to the linear regime. It would certainly be interesting to 
pursue the calculations using numerical simulations as one could 
expect that the holes produced in the plasma would widen with 
the growing magnetic pressure inside and that the density excesses 
become spikes. Thus one might reasonably expect the generation 
of many local small magneto-sonic shocks that are absorbed by 
the main ultra-relativistic shocks, with the reconversion of some 
amount of cosmic-ray energy into thermal energy in the precursor. 

The previous discussion has focused on the role played by the 
first generation of cosmic rays, since higher energy cosmic rays 



P..o 



(42) 



where k,^ denotes the same wavenumber evaluated for the whole 
cosmic-ray distribution, and p,o denotes the minimum momentum 
of the cosmic ray distribution at the shock front. Consequently, if 
7j- oc ^, with a = 2/3 or 1/2 for the two branches of the instability 
obtained above: 

7,4... = -^ rv.o4, , (43) 

\P:qI 

with b = I for large scale turbulence, b = 2 for short scale tur- 
bulence. Here as well, 7,0 should be understood as the growth 
rate derived previously for the whole cosmic-ray distribution. This 
product grows with increasing p,/p,.o provided s < (b + a) I a, i.e. 
s <%{b = 2,a = 2/3), ox s < 5 {b = 2, a = 1/2), or s < 5/2 
{b = 1, a = 2/3) or finally s < 2i {b = \, a = 1/2). These inequal- 
ities are likely to be fulfilled. Finally, the product A;,/mhd scales as 
(p,/p,,o)' ' hence decreases if s > 1. This also means the typical 
scale of the instability increases with increasing p,/p,,o. 

Of course, the above implicitly ignores the spatial dependence 
of 6, but it suggests that amplification through the streaming of the 
highest energy part of the cosmic ray distribution will seed much 
more efficiently the magnetic field amplification. 

3.1.5 Non-existence of an incompressible instability 

Since the above compressive instabilities appear to saturate at a 
level too small (SB/B ~ 1) to allow Fermi acceleration over a broad 
range of energies, the possible existence of non-compressive insta- 
bility becomes crucial. One can search for such modes by selecting 
the wave vectors accordingly, in which case the reduced system 
becomes: 



D^u. - I ^ + 



k,dz 



u, k,dy \ k, 

kl /SshFsh k] j ' k. 



k,d. 



Du,- 



k.d] 



a. ~, 

—D-u. 

k, ~ 











(44) 



It is easy to verify that this system possesses no unstable solution 
but only damped solutions. Hence incompressible unstable modes 
do not exist in this MHD description. 



4 CONCLUSIONS 

In this study, we have examined the amplification of a pre-existing 
magnetic field upstream of an ultra-relativistic shock wave. We 
have assumed that the magnetic field is fully perpendicular to the 
shock normal in the shock front frame, as generally expected in the 
ultra-relativistic limit Fsh » 1. In the shock frame, the cosmic rays 
do not induce any current at zeroth order, only a net charge den- 
sity, which is partly screened by the inflowing plasma. This charge 
distribution then triggers an instability on very short spatial scales. 
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with a growth rate increasing with the wavenumber. Cosmic rays do 
not respond to this amplification as it takes on scales much shorter 
than the typical Larmor radius. 

It is important to stress that these instabilities are of a differ- 
ent nature than the resonant or non-resonant instability studied by 
Bell (2004) in the case of a non-relativistic parallel shock wave, 
most notably because they are essentially compressive. Since den- 
sity depletions are naturally limited, our linear analysis indicates 
that these instabilities saturate at a moderate level of amplifica- 
tion, 5S/So ~ 1 . Furthermore, there is no unstable incompressible 
mode within the present MHD approximation which could provide 
a higher level of amplification. Therefore, we conclude that, within 
the framework of ideal MHD, instabilities at ultra-relativistic mag- 
netized shock waves appear limited in their efficiency. 

Such instabilities cannot therefore account for the degree of 
amplification that has been inferred from the modeling of the af- 
terglow radiation of gamma-ray bursts. We have also argued in 
Section 12.21 that the ratio SB/Bq sets the dynamic range of ener- 
gies over which powerlaw spectra can be produced through Fermi 
acceleration. The present instabilities thus cannot allow success- 
ful Fermi acceleration. Nevertheless, it would certainly be useful 
to pursue the present investigation with dedicated numerical sim- 
lations, which would allow one to go beyond the present linear 
approximation and study whether non-linear effects can push the 
small scale magnetic field to higher values. 

If the standard interpretation of gamma-ray burst afterglows as 
the synchrotron light of electrons accelerated at the forward shock 
(but see Uhm & Beloborodov 2006, Genet, Daigne & Mochkovitch 
2007 for recent alternatives involving the reverse shock) holds, one 
is led to conclude that some other instability is able to amplify 
the upstream magnetic field on short spatial scales, or that some 
other form of acceleration mechanism is operating (see for instance 
Hoshino 2008 for an alternative involving radiative pressure ef- 
fects). The present work indicates that such cosmic ray induced in- 
stabilities would involve non-MHD effects. The short spatial scales 
Ac <K TL/Fsh required are in conflict with the usual synchrotron res- 
onance between cosmic rays and MHD waves; this remark is actu- 
ally one of the motivation of the present work in a full MHD frame- 
work. However, one cannot rule out that other types of resonance 
could become relevant and yet involve incompressible modes, such 
as Alfven waves. One promising instability, currently under study, 
involves a Cerenkov resonance (co = k^c) between plasma waves 
and the cosmic ray beam, generating a modified two stream insta- 
bility. The growth rates appear quite promising and the details will 
be given in a forthcoming article. 



APPENDIX A: BACKGROUND AND PERTURBED 
QUANTITIES FOR PERPENDICULAR SHOCK WAVES 

Al Stationary regime 

In the shock front frame, the electric field, current and charge den- 
sity read: 

1 U-By UyBy 

E - — u X B = ^— ex ^ez , 

j = — <9vfivez , 
An 



Ppi 



-Pci + r-5v£v 
An 



(Al) 



The four velocity of the inflowing plasma is written (yc, uc) (hence 
u = y/3 is dimensionless). In this equation, y represents the total 



Lorentz factor of the upstream plasma; it will be shown to coincide 
nearly with Fsh in the following. 

Note that in the shock front, there is no current at zeroth order 
associated with the cosmic rays, only a net charge density p„ which 
may be partially screened by an induced non-zero plasma charge 
density Ppi. Faraday's law implies in this stationary regime: 



V X E = ^ <9.t 



UyBy 

7 



0. 



Hence the Gauss equation can be rewritten as: 



Ppi 



+ \OxU, a.jM.v 



(A2) 



(A3) 



Any \ ' Kj- 
Energy momentum conservation for the upstream fluid im^ 



plies: 



9 1 

p„c^u ■ Vu = pp,E + -J X B , 



(A4) 



where p„ is the proper mass density of the upstream plasma. De- 
noting Fm = -puC^Mv the mass flux (times c) through the shock 
front: 



FmdxUx 



FmdxUz 



OjM; OxUy 



Any \ 
47r ' 



\t - -; \OxUz 0.v".v 



UyBy 

7 



Any \ 
Since: y = {\ + ul + ujy^, 

0x7 = —— '■ — 

7 

therefore ( IA2t implies: 



y I y^ — u\ UxU, 

dyBy = -By T—SyUy T^d^U- 

Ux \ y^ 7 



Using l lA7t . the system JA5t can be transformed into: 



(A5) 



(A6) 



(A7) 



Bl 



Any^Ux 
Blux 



dyUx 



Any^ [ Any^ 



Pcr- 



-Pc- 



U^By 



7 

UxBy 

7 



(A8) 



This sytem can be solved by successive approximations. For 
instance the ratio: 



B. 



Any^UxF„ 



Fr/31 « 1 



(A9) 



since \Ux\ ~ r^i^;/^^ = B,!,,/ ^jAnp^c- represents the Alfven velocity 
(as measured in the upstream rest frame). Hence: 



dyUx - Pci 



U-By 

7Fm 



(AlO) 



Consider now the ratio of the first term on the Ihs of the second 
equation in JASt to the term of the rhs of the same equation: 



dxU> 



Any- 



Pcr- 



llxBy 






(All) 
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This ratio is much smaller than unity (since u- <c Fji,). Then: 



drU- 



Pci 



f + -2. 



^^(x) 






p„dx . 



(A 12) 



In the last equation, x, corresponds to the maximal distance to 
which cosmic rays can stream ahead of the shock front, as mea- 
sured in the shock front frame. We also assume that j(„ 6,, and p 
vary on much longer scales than u,, and that \u.\ « |mv|. This hier- 
archy can be verified by taking the ratios of the equations in ( |A5t . 
Since BlliAny) » P\F^, one can simplify the previous ex- 
pression and obtain the following order of magnitude: 



\Uz\ 



4,- 



rL, TshPuC^ 



(A13) 



To obtain this result, we have approximated the integral over p„ 
as 4rPcr, 4r representing the length scale of the distribution in the 
shock frame, and written p„ - e^J{p,c); e^r denotes the cosmic ray 
energy density while p, (r^.) represents their typical momentum 
(Larmor radius). Finally, ecr/(rshPuC^) ~ ^A^a where ^cr is the 
fraction of shock internal energy carried away by the accelerated 
population. In a straightforward way, one obtains: 



'L. 



-0(1), 



so that: 

fci <; 1 ■ 

\lh\ ^ 

This justifies the previous approximations. 



(A14) 



(A 15) 



A2 Time dependent quantities 

A2.1 Perturbed electric field 
The perturbed Ohm law leads to: 
1 



SE 



- BySu- + u-SBy + 
yl ~ 

+ — (rshj8sh"c(5M.,Bv - Su,ulByj ex -I- 

- {-SB^u- + SB^u^) ey + 
7 

- TshAhfJSv - (5k„Sv + 
yL 



(A16) 



We keep track of the terms in u- in these equations in order to 
properly evaluate the terms which involve the derivative of u- with 
respect to x. This latter quantity is of the same order than other 
background quantities and cannot be neglected. However, terms in 
u- will be neglected in the final equations of evolution of the various 
background quantities since u- <K m^-. For similar reasons, deriva- 
tives of background quantities (other than u-) such as Ppj, By and 
Uj (hence y) can be neglected. Thus, setting y ^ Fsh, /3 - ySsh we 



obtain: 



6E 



{:^Bydu^ + -—U-6B,. + -^u-BySu^ 

■* ''^ sh 

:)ex ■ 



H Yu':BrSii-]e^ + 

^ sh 



A sh 



(/3,h5Sy - -^B,li,^Su,u, - --S,<5M,)e, . (A17) 



A2.2 Perturbed equations 

The perturbed equations for the perturbed components of the mag- 
netic field, defined in units of By as b = (5B/6v can be written: 



bb^ = -^dySux + uA-d-b^ + -^dySu^ 
Dby = -d-6u- t-3j(5mv - -^u'.du- + 

+u-\—d-by —d-diiy —dySu- + 

"^ " ' Fsh ^ Fsh 

Db^ = dySU^ + U'.by + U-idyby + dyby H —dySUy + 

~\ Fsh 

-^dydu]. (A18) 

1 sh 

The continuity equation can be expressed as: 



^ Sw 1 

D — = -V -Su d,6y - -(dySuy + dySuy + d-6u,) . (A19) 

w c . . , , 



We have noted w = p„c^ the proper enthalpy density of the up- 
stream plasma. In the perturbation regime the contribution of the 
temporal derivative of the fluctuating part of the Lorentz factor is 
negligible. Again, while deriving Eg. d A 19b the x derivatives of ze- 
roth order quantities have been discarded. 

The equation of motion in the shock rest frame accounting for 
the displacement current are: 



d2 o 1 ^u 

wDu=PpiE-V— +B-V xB. (A20) 

^ 8/r An Anc dt 



Using Eq. dAlVI ) and developping the Ea. JA20t to the first order we 
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finally get after some manipulations: 

bSu, = 0iTl^ [dyb, - d,b,) - kSu, + 

+131^1 (pJ-dA - :^d,6u\ - p%— + 



+ M, 



A* 



A sh 

9 Sw / \ 

DSuy = (Jsh^shKb- - P:dy h M, l/fi'v - d-dUy) , 

■ ■ vv ~ ^ ~ ■ / 

[l+Pl] D6u, = Pl [dyb, - d-M) - Ahrsh/^ ( 1 - /3i) foy + 
-<(l-/?i)&,-^/3i5>,+ 

1 sh 



6w 



6w 



-Ptd,— +PaJ'AK— + 
w w 



+11- 



ySsh 



A sh 

+P.bT,^l [dyb, - d,A) + 

+PshP'A -d,6Ux - AhS.v^M.v + 



+l3iT,A-d,by-P.i.d,by 



(A21) 



We used the shorthand notation Pi = c^/c^ with fj the speed of 
sound, and: 



Pi^ 



4/rr> 



(A22) 



which defines the Alfven velocity squared (as measured upstream). 
The quantity k, which carries the dimension of a wavenumber, is 
defined as follows: 



1 Pcrfiy 

l+Pl TshW 



(A23) 



The plasma charge density ppi is indeed related to the cosmic-ray 
charge density Pd through the Maxwell equation: 



Ppl = -Pcr + 



By , 



47rnt 



(A24) 



hence: 



Ppi By 

TshW 



(A25) 



evolution of the perturbed quantities reads 
1 



Db, 
Dby 

Db, 

DSUy 
DSUy 

bdu- 



w 



dySuy: + u- -d.b^ + ——dySuy 






-d-Su- —dySuv —kSu, + 

+uA-d-by - -^d-Suy - -^dySuA , 
dySuz + Kb J + u-{—d-b- + -^dySiiy 1 , 
-k6u, - w. /cfc,, +j8sii — 6uy + d-6uA , 

Psh^shKb- + U, [Kby — d-Sllyj , 

-pshTshKby + —rSuy + pshTshK — + Uz -^kSu-^ - dz6u- 1 , 



-dj6uj - dySuy - dz6u. 



(A26) 



To analyze the instabilities, it is convenient to make the fol- 
lowing reduction: D = kyP^h^shb, Suj = pstX^hiij', bj and — un- 
changed. The ordering for doing simplifications is consistent with 
the expected physics: we assume ky,k,,k, <c ky, u- <K P^b^sh and 
Du, = -TshK/ky at the same order than previous terms. The sim- 
plified system reads: 



Db, 


~ 





Db,. 


^ 


dz k, , 
~"^-P.J-^u,-Pl. 


Db- 


^ 




bih 


^ 


1 k. 


blly 


^ 


±'-ib. 

P^k, - 


Du- 


^ 


1 k, 1 k,6w 

Psh h " ySsh kj, w 


-Sw 




/a,- a,, a,. 


D — 

w 


— 


-[V'T'T 



Uz 5.V . 

«- 

Psh^A kx 



(A27) 
(A28) 

(A29) 

(A30) 

(A31) 

(A32) 

(A33) 



Thus one gets the intermediate system: 



D'ff, 



D^tL = 



fc, k,d, u, ktdx\ k, - Sw 

kj ySshKj Msh kj. I Pshkx W 



k,dx 



Du- - 



k,dl 



i^z -2 

u. — -—D u. 



P,i,kl '- p,^kl ' k. 



(A34) 



In the limit ySA « 1 and /3s « 1, the system governing the 
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